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We demonstrate that a reduction in the noise-to-signal ratio may be obtained for hadrons at non- 
■ zero momenta whilst maintaining a good overlap with the ground state through a generalisation 

of Gaussian/Wuppertal smearing. The use of an anisotropic smearing wavefunction is motivated 
_C by the physical picture of a boosted hadron. 
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1. Introduction 



Figure 1: Effective energies for pions at different values of q 



Lattice QCD is a useful tool for the extraction of quantities such as masses, energies and 
matrix elements. These may be extracted for a given hadron from the Euclidean time dependence 
of correlation functions for states with appropriate quantum numbers. The ground state may be 
isolated through taking the large source-sink separation limit due the time dependent exponential 
suppression of excited states. However the statistical errors of the correlation functions also grow 
with increasing time separations (with the exception of the pion at zero momentum), and so it can 
prove difficult to find a 'window' where both statistical and systematic errors are under control. 
This becomes increasingly more difficult once states with non-zero momenta are considered, as is 
the case for many desirable quantities such as form factors. A clear example of the deterioration 
of the signal for increasing momentum q 2 is given for the pion in figure 1. Throughout this paper 
we will use the notation p = q ^ when referring to the momenta earned by the hadron, where 
q = (<7i,<72,<?3) is a vector of integers. For the case of the pion at zero momentum, the noise-to- 
signal ratio is expected to be constant as a function of time whereas for non-zero-momentum it is 
expected to be [1, 2] 

Rns(xo) « e ( V»4+p 2 -»a*o ? foTXQ x _ (L1) 

This behaviour is demonstrated for our data in figures 1 and 3 for several values of q 2 . The asymp- 
totic trend sets in quite early and is consistent with noise dominated by a zero-momentum two pion 
state. 

The usual technique used to improve the ground state overlap with the interpolating source 
fields, thus improving the 'window', is through smearing. The present study is devoted to opti- 
mising the smearing based upon the hadrons kinematics. In the following sections we detail how 
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Figure 2: The yz-section of a spherical wavefunction K x = K y = K z = 2.9 (left panel) and an anisotropic 
wavefunction K x = K y = 2.9, K z = 2.9/16 (right panel). 

a generalisation of Gaussian/Wuppertal smearing that allows for anisotropic quark smearing (in- 
terpolating a boosted hadron) can lead to a significant improvement of the noise-to-signal ratio for 
states of non-zero momentum. The calculations were performed using Nf = 2 non-perturbative 
&(d) improved Wilson fermions on configurations generated as part of the CLS initiative. This 
initial exploratory calculation utilises 168 configurations, each with four evenly spaced source po- 
sitions, on a lattice of 32 3 x 64 points, with a lattice spacing of 0.063(2) fm and a pion mass of 
approximately 450 MeV; further details on the simulation parameters may be found in [3, 4, 5]. 

2. Implementation of anisotropic smearing 

Gaussian/Wuppertal smearing is defined by [6, 7] 

H k (x,x') = VkWS^ + Ufa-ktf^t , (2.1) 
where H k is the hopping operator in direction k and the smeared fermion field xpW is defined as 

¥ W(x) = {\ + KH)"y(x), (2.2) 

where 

H = H X +H y +H z . (2.3) 

In our generalization, we promote K and H to spatial vectors K and H with components x, y and z 
so that 

Y {K) (x) = (l + r-H)*y(x) • (2.4) 

It is then possible to interpolate a boosted hadron through tuning the components in K so that a 
Gaussian-like wavefunction is squeezed into something more point-like in the boosted direction. 
We also note that this tuning between point-like and gaussian-like enables a balance of parameters 
that result in, in the first case a better noise-to-signal ratio and in the latter a better overlap and thus 
plateau [8]. Figure 2 shows the jz-cross sections corresponding to K x = K y = K z = 2.9 (spherical - 
left panel) as well as K x = K y = 2.9, K z = 2.9/16 (anisotropic - right panel). We use fQ = 2.9 as a 
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Figure 3: Noise-to-signal ratio for a pion for different q 2 plotted against the expected asymptotic behaviour 
for anisotropic smearing. The grey lines show the corresponding behaviour for spherical smearing (the 
corresponding spherical data points are omitted). 



basis and fix n = 140 as these values maximised the length of the effective mass plateaus for the 
nucleon and pion at q 2 = and produce a wavefunction with a radius of approximately 0.5 fm [4]. 
The anisotropic smearing is implemented so that only directions with momenta are squeezed. In 
addition we smear the gauge links in the hopping operator by applying one level of hypercubic 
smearing, using the parameters specified in the definition of "HYP2" links in [9]. 

3. Pion and nucleon effective-mass study 

We consider two-point correlation functions of the form eq. (3.1), where the operator O inter- 
polates mesons or baryons 

C(x) = (o(^ K \x),W {K \x))O\^ K X0),W {K) m) ■ (3-D 

Through the use of the anisotropic smearing we aim to reduce the contribution to the noise from the 
zero-momentum component of the four-point correlation function occurring in the variance of the 
two-point correlation function. Figure 4 demonstrates the effect of anisotropic smearing for four 
values of q 2 . For each direction that the squeezing is applied to, the boost factor y, defined by 

Y= X , (3-2) 

m H 

is indicated. 
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Figure 4: Effective energies for the pion for different values of q 2 using spherical (blue) and anisotropic 
(red) smearing. 



The results demonstrate that anisotropic smearing offers a clear improvement for non-zero mo- 
menta, which becomes more pronounced for larger boosting factors. This is further demonstrated 
in figure 3. We investigated the effects of adjusting the anisotropy in different directions by multi- 
plying the appropriate components of K with a factor 1 /N, where N = 2, 4, 8, 16 and 32. We found 
that by Af = 16 the improvement in the statistical noise had saturated and any further improvement 
was accompanied by a slightly worsened overlap with the ground state. We explicitly checked that 
squeezing in directions orthogonal to those with momenta provided no gain. The same techniques 
were applied to the nucleon (figure 5), however, due to the worse initial signal and larger mass, the 
nucleon only attains relatively small boosting factors before the signal is lost in noise and so the 
effect is less pronounced. 

Figure 6 compares the extracted results for the spherical and anisotropic smearing to the con- 
tinuum dispersion relation for both the pion and nucleon, further demonstrating the potential sta- 
tistical gain. The extracted values were obtained via a % 2 per d.o.f. minimisation of simple plateau 
fits to the data in figures 4 and 5, and the same fitranges are used for both spherical and anisotropic 
smearing. 
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Figure 5: Effective energies for the nucleon for different values of q using spherical (blue) and anisotropic 
(red) smearing. 




Figure 6: The dispersion relation plotted against the extracted effective energies for spherical (blue) and 
anisotropic (red) smearing for the pion (left panel) and nucleon (right panel). The solid line denotes the 
continuum dispersion relation. 
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4. Conclusion and outlook 

We have implemented a generalisation of Gaussian/Wuppertal smearing which allows the 
smearing to be optimised to account for a hadron's kinematics, which for a boosted hadron cor- 
respond to an anisotropic wavefunction. Our data demonstrates that this provides a significant 
improvement in the noise-to-signal ratio for heavily boosted hadrons whilst maintaining a good 
overlap with the ground state. The results are more pronounced for the pion due to its large boost- 
ing factor at small values of the lattice momenta. Our results show that the improvement is more 
dependent on y than on the mass of the hadron, and as such we expect this prescription to be more 
effective as the chiral limit is approached. Important applications for this method could include 
the determination of mesonic and baryonic form factors, such as for B — > %lv or the nucleon form 
factors Ge and Gm, for which we presented preliminary results at this conference [10]. A prelimi- 
nary study of the nucleon three -point function used in the calculations of the nucleon form factors 
show that a similar gain could be achieved. We are currently extending this study. In principle, 
anisotropic smearing requires a different smearing for each choice of momentum and so the largest 
gain would likely be achieved though the use of all-to-all [11, 12] stochastic propagators to reduce 
the number of inversions of the Dirac operator and thus the computational cost. 

A published account of this work may be found in [13]. We note that a similar approach was 
applied to the nucleon channel in [14], which is consistent with our findings. 
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